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ABSTRACT 

A theoretical light-curve model of the 1985 outburst of RS Ophiuchi is presented based on a ther- 
monuclear runaway (TNR) model. The system consists of a very massive white dwarf (WD) with an 
accretion disk (ACDK) and a red giant (RG). The early phase of the V light curve is well reproduced 
only by the bloated WD photosphere of the TNR model on a 1.35 ±0.01 Mq WD while the later phase is 
dominated both by the irradiated ACDK and by the irradiated RG underfilling the inner critical Roche 
lobe. The UV light curve is also well reproduced by the same model with the distance of 0.6 kpc to RS 
Oph. The envelope mass at the optical peak is estimated to be 2 x 1O~ 6 M , indicating a rather high 
mass accretion rate of 1.2 x 10 _7 M Q yr _1 between the 1967 and the 1985 outbursts. About 90% of the 
envelope mass is blown off in the outburst wind while the residual 10% (2 x 10~ 7 M Q ) has been left and 
added to the helium layer of the WD. The net increasing rate of the WD mass is 1.2 x 1O _8 M0 yr _1 . 
Thus, RS Oph is certainly a strong candidate for a Type la supernova progenitor. 

Subject headings: binaries: close — binaries: symbiotic — novae — stars: individual (RS Oph) — 
supernovae 



1. INTRODUCTION 
RS Oph is one of the well observed recurrent n ovae and 



characterized by a lo ng orbital period of 460 days ( Dobrzy 



cka & Kenyon 1994), a relatively short recurrence period 
of ~ 10 — 20 yrs, and a companion (M-giant) star under- 
filling the R oche lobe but losing its mass by massive stellar 
winds (e.g., Dobrzycka et al. 1996). Historically, RS Oph 
underwent five outbursts, in 1898, 1933, 1958, 1967, and 
1985, with t he light curves very similar each other (e.g. 



Rooino 1987). The latent 1085 outburot hao boon oboorved 



at all wave lengths from radio to X-rays (e . g ., papers in t h e 



1985 RS Oph conference proceedings ed. by |Bode 1987| ) 
Although there had been inten se debates on the mecha- 



nism of RS Oph outbursts (e.g., Livio et al. 1986b| |Web- 
bink et al. 1987| ), various observational aspects favor ther- 



monuclear runaway (TNR) models on a very massive white 
dwarf (WD) (e.g., Anupama fc Mikolajewska 1999 for a 



recent summary). Rapid decline rates of the light curves 
indicate a very massive WD close to the Chandrasekhar 
limit. Kato (1991) has first calculated RS Oph light curves 
for the WD masses of 1.33, 1.35, 1.36 and 1.37 M Q and 
found that the light curve of the 1.36M Q model is in better 
agreement with the observational light curve of RS Oph 
than the other lower mass models. 

However, structures of bloated env elopes on white 
dwar f s have bee n drastically changed ( Kato fc Hachisv] 
1994 ; Kato 19991) after the advent of the new opacity (e.g. 



Iglesias fc Rogers 1996 ). Thus, we have been recalculating 



light curves of the recurrent novae and have again deter 
mined the mass, com position, and distance ( Hachisu & 
Kato 19991 for T CrB; lHachisu et al. 2000ai 2000b for U 
Sco; |Hachisu fc Kato 2000| for V394 CrA). Moreover, we 



have added new parts to Kato's (1991) simple model in 



order to follow later stages of the light curves, i.e., the 
irradiation effects of the red giant (RG) companion star 
and the accretion disk (ACDK). As a result, much more 
reliable quantities are derived from the light curve fittings 
than Kato's ones. 

Type la supernovae (SNe la) are one of the most lumi- 
nous explosive events of stars. Recently, SNe la have been 
used as good distance indicators which provide a promis- 
ing tool for determining cosmological parameters because 
of their a l most uniform maximum luminosities ( Riess et 
al. 19981; |Pcrlmutter et al. 1999| ). These both groups 
derived the maximum luminosities (£ max ) of SNe la com- 
pletely empirically from the shape of the light curve (LCS) 
of nearby SNe la, and assumed that the same L max -LCS 
relation holds for high red-shift SNe la. To be sure of any 
systematic biases, the physics of SNe la must be under- 
stood completely. By far, one of the greatest problems 
facing SN l a theorists is the lack of a real progenitor (e.g., 
Livio 1999 for a recent review). Finding a reliable pro- 
genitor is urgently required in SN la research. Recurrent 
novae are probably the best candidate for this target (e.g., 
Starrfield, Sparks, & Truran 1985; Hachisu et al. 1999b; 
Hachisu, Kato, & Nomoto 1999a). 

In §2, we present our methods for calculating theoretical 
light curves based on the optically thick wind theory de- 
veloped by Kato & Hachisu (1994). Fitting results of the 
modeled light curves to the observations are shown in §3. 
Our results indicate that RS Oph is a strong candidate for 
Type la supernovae. Discussion follows in §4 especially for 
the distance to RS Oph because we obtain a rather short 
distance of 0.6 kpc. 

2. THEORETICAL LIGHT CURVES 
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Light Curves of RS Ophiuchi Outbursts 



In the TNR model, it has been established that the WD 
photosphere expands to a red giant size and then gradually 
shrinks to the original size of the WD in quiescent phase 
(e.g., - 0.0039 R Q for M W d = 1.35 M Q ) with the bolo- 
metric luminosity being kept near the Eddington limit. 
The visual luminosity reaches its maximum at the max- 
imum expansion of the photosphere and then gradually 
decays to the level in quiescent phase. Accordingly, the 
main emitting region of the WD photosphere moves from 
visual into soft X-rays through ultraviolet. Optically thick 
winds, which are blowing from the WD during the out- 
bursts, play a key role in determining the nova duration 
because a large part of the envelope mass is blown in the 
outburst wind. The development of the WD photosphere 
can be followed by a unique sequence of steady-state, op- 
tically thick, wind solutions as shown by Kato & Hachisu 
(1994). 

We have calculated such sequences of Mwd = 1-3, 1.35, 
1.36, 1.37, and 1.377M by assuming the hydrogen con- 
tent X = 0.70, 0.50, 0.35, the helium content Y = 0.28, 
0.48, 0.63, respectively, and the metallicity Z — 0.02 (so- 
lar metallicity) and then, obtained the optical li ght curves. 
Here, we hav e used the updated OPAL opacity ( Iglesias & 
Rogci ^l996 ) . The envelope mass is decreasing due to wind 
mass loss (M w ind) and hydrogen shell burning (M nuc ), i-e., 



dt 



-AM = M„. 



(1) 



where M acc is the mass accretion rate of the WD and as- 
sumed to be M acc — 1 x 1O _7 M yr^ 1 during the outburst. 
The photospheric temperature T p h, the photospheric ra- 
dius i? p h, the photospheric velocity t> p h, the wind mass 

loss rate Af w ind> and the nuclear burning rate M nuc are 
unique functions of the envelope mass AM. Integrating 
equation ([!]), we can follow the development of the enve- 
lope mass AM and obtain various physical quantities of 
the WD envelope. The mass lost by the wind, AM w m, 
and the mass added to the helium layer of the WD, AMh c , 
are calculated from 

AM wind = J M wind dt, AA/hc = J M nuc dt. (2) 

Assuming a black-body photosphere, we have calculated 
the ^-magnitude of the WD photosphere with a response 
function given by Allen (1973). For simplicity, a limb- 
darkening effect is neglected. By fitting with the observa- 
tional point s in the early 4 days of the superposed V light 
curve (e.g., Rosino 1987 ) and the 1985 V light curve, we 
have determined the WD mass of Mwd = 1-35 ± 0.01M Q 
and the apparent distance modulus of too = 11.09 as 
shown in Figure |l|. The early 4 days ^-magnitude is de- 
termined mainly by the bloated WD photosphere, and de- 
pends sensitively on the WD mass, but depends hardly on 
the hydrogen content, X, the binary parameters such as 
the mass ratio, or the accretion disk shape, so that the 
WD mass determination itself is rather robust (see also 



Kato 1999| ). 

To reproduce the V light curve in the late stage (t ~ 4 — 
100 days after the optical maximum), we have to include 
the contributions of the irradiated RG and the irradiated 
ACDK. As suggested by Dobrzycka et al. (1996), we have 



assumed that the RG lies well within the inner critical 
Roche lobe, i.e., 



Rkg^jR; ( 7 <i), 



(3) 



where Ri, is the effective radius of the inner critical Roche 



lobe for the RG component (e.g., Egglcton 1983) and 7 is 
a numerical factor. Here, we assume 50% efficiency of the 
irradiation (?7rg = 0.5). The nonirradiated photospheric 
temperature of the RG is a parameter for fitting, and has 
been determined to be T p h,RG = 3100 K for 7 = 0.4. 

The orbit of the companion star is assumed to be circu- 
lar. The ephemeris for t he inferior conjunction of the giant 
is 2,444,999.9+460 x E ( pobrzycka fc Kcnyon 1994| ). The 
light curves are calculated for six cases of the companion 
mass, i.e., M RG = 0.5, 0.6, 0.7, 0.8, 1.0, and 1.2M . Since 
we obtain similar light curves for all of these six masses, 
we show here only the results for Mrq = O.7M . In this 
case, the separation is a — 318.5i? , the effective radii of 
the inner critical Roche lobes for the WD component and 
the RG component are R\ = 139. LR S and Rt, = 103.1i? o . 
If 7 = 0.4, then we have R 2 = OAR* ~ 4Oi? . 

We have also included the luminosity coming from the 
ACDK irradiated by the WD photosphere when the accre- 
tion disk exists during the outburst. Here, we assume 50% 
efficiency of the ACDK irradiation (?7dk = 0.5). The vis- 
cous heating is neglected because it is rather smaller than 
that of the irradiation effects. The temperature of the un- 
heated surface of the ACDK including the rim is assumed 
to be Tdisk = 2000 K. We have checked two other cases of 
Tdisk = 1000 and K and have found no significant dif- 
ferences in the light curves. The shape of the ACDK is 
assumed to be axisymmetric and approximated by 



-Rdisk = aR\, 



and 



h = (3R, 



disk 



-Rdisk 



(4) 



(5) 



where i?disk is the outer edge of the accretion disk, h the 
height of the surface from the equatorial plane, w the dis- 
tance on the equatorial plane from the center of the WD, 
and a and (3 are both numerical factors which should be 
determined by fitting. The power of v is assumed to be 
v = 9/8 from the standard disk model. We have checked 
the dependency of the light curves on the parameter v 
by changing from v = 9/8 to v = 2 but cannot found 
any significant differences when the disk rim is as small as 
/3 = 0.01—0.05. 

3. RESULTS 

The early 4 days V light curve can be reproduced with 
the WD mass of Mwd = 1-35 ± O.O1M as already men- 
tioned in the previous section. The late phase V light curve 
(t ~ 4 — 100 days after maximum) indicates strong irradi- 
ations both of the ACDK and of the RG. If we assume a 
lobe-filling RG companion (7 = 1), then the irradiation of 
the RG is too luminous to be compatible with the obser- 
vational light curves. On the other hand, a smaller size 
of the RG such as 7 = 0.4 (~ 4Oi? ) gives a reasonable 
fit with the observational one as shown in Figure 0. Here, 



we assume th e inclination angle of i — 30° (Dobrzycka & 
Kenyon 1994}). 
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To fit the late phase light curve, we have finally adopted 
the disk parameters of a — 0.01 and (3 = 0.01 as shown 
in Figure |[ Here, the thick and thin solid lines denote 
the cases of (a = 0.01, (3 = 0.01) and (a = 0.1, = 0.01), 
respe ctively. Something between these two can roughly re- 



produce the light curve of the 1985 outburst but the case 
of (a = 0.01, j3 = 0.01) is much better to reproduce the 
superposed light curve given by Rosino (1987). 

We have also calculated UV light curves with a response 
function of 91lA— 3250A to fit the UV data flBnijder;" 



1987a). Figure || shows three cases of UV light curves for 
X = 0.35, 0.50, and 0.70, with the parameters of a = 0.01 
and j3 = 0.01. The calculated UV light curves can re- 
produce well the UV observations for X — 0.70 and the 
distance to RS Oph of d = 0.57 kpc. Thus, we have ob- 
tained the absorption of A v = 11.09 - 5 log(570/10) = 2.3 
for d = 0.57 kpc. This absorption of Ay — 2.3 is consistent 
with the color excess of E(B - V) = 0.73 flSnijdcrs 1987b] ) 
because of Ay = 3.1 E(B — V), although this distance 
of ~ 0.6 kpc is not consistent with another estimation of 
1.6 kpc from the hydrogen column density (Hjellming et 
al. 19|(]). The UV light curves depend hardly on the disk 
parameters of a, f3, or v when a ;> 0.01 mainly because 
the UV light is coming from the innermost part of the 
ACDK, but depend on the irradiation efficiency of ?7dk- 
We have examined that the distance becomes 0.70 kpc for 
rjDK = 1-0 or becomes 0.50 kpc for ?7dk = 0.25. 

The envelope mass at the optical maximum is esti 



matec to be AA7 = 2.1 x lO^Ai^, which is indicating 
a mass accretion rate ot O x iu Mq yr 1 during the 
quiescent phase between the 1968 and the 1985 outbursts. 
About 90% of the envelope mass (1.9 x 1O~ 6 M ) has been 
blown off in the optically thick wind and the residual 10% 
(0.2 x 10~ 6 M Q ) has been left and added to the helium 
layer of the WD. The residual mass itself depends on both 
the hydrogen content X and the WD mass. It is roughly 
ranging from 7% (X = 0.70) to U%(X = 0.50) including 
ambiguities of the WD mass (M WD = 1.35 ± O.O1M ). 
Therefore, the net mass increasing rate of the WD is 
1.2 x 10~ 8 M Q yr" 1 , which meets the condition for Type 
la supernova explosion if the WD core consists of carbon 



and cjxygen. Thus, we conclude that RS Uph is an im- 
mediate progenitor of Type la supernova if the donor is 
massive enough to supply fuel to the WD until the WD 
reaches 1.378M Q for explosion (Nomoto, Thielemann, & 
Yokoi 1984). 

4. DISCUSSION 

Assuming a tilting accretion disk around the WD, 
Hachisu & Kato (1999) have reproduced theoretically the 
second peak of T CrB outbursts. Such a radiation induced, 
tilting disk instability sets in if the condition 



Ma 



3 x 10- 8 M Q yr 



/ -Rdisk 

- 1 55 I R 



1/2 




i?WD 



0.004 R 



2 x 10 38 erg s -1 

(6) 



M\VD 
1.35 Mr. 



is satisfied ( Southwell et al. 1997 ). Our estimated ac- 
cretion rate of the WD in RS Oph is about M acc ~ 
1.2 x 1O~ 7 M yr" 1 , which does not meet the above con- 
dition, so that we do not expect the radiation induced 



instability in RS Oph system. Thus we can explain the 
reason why a second peak as seen in T CrB system does 
not appear in RS Oph system. 

Very soft X-rays were observed 251 days after the opti- 
cal maxi mum in the 1985 outburst of RS Oph (Mason et 
al. 1987j ). Mason et al. estimated a black-body tempera- 
ture of 3.5 x 10 5 K, a total energy flux of L x ~ 1 x 10 37 
erg s _1 (d/1.6 kpc) 2 for a hydrogen column density of 
Nh = 3 x 10 21 cm -1 , and concluded that these soft X- 
rays come from a WD photosphere with a steady hy- 
drogen shell-burning. However, steady hydrogen shell- 
burning has stopped 122 days after the optical maximum 
in our model of X = 0.70 and M WD = 1.35M (Fig. |), 
thus indicating an accretion luminosity instead of hydro- 
gen shell-burning. Using the black-body temperature of 
T ~ 3.5 x 10 5 K and the WD radius of 0.0039 R Q , we esti- 
mate a total luminosity of Lx = ATiR^ NT) aT A ~ 2OOL0. 
This low luminosity corresponds to the accretion lumi- 
nosity of = 2L x Rwd/GM wd ~ 0.5 x 10~ 7 M Q 
yr -1 , which is roughly consistent with our estimated value 
of M acc - 1.2 x 10~ 7 M Q yr -1 . Then, a rather short 
distance of 0.45 kpc to RS Oph can be derived from 
8 x 10 35 - 1 x 10 37 (d/1.6 kpc) 2 . Thus, our short distance 
of 0.6 kpc is consistent with the soft X-ray observation. 
Moreover, the recent optical and ultraviolet observations 
also s uggest 100 — 600 L Pl fo r a hot component of RS Oph 
(e.g., Dobrzycka et al. 1996 ). 

X-r ays were also observed at quiescent phase ( |Oric 
1993), but the flux is too low to be compatible with the 



TNR model. One possible ex planation is an abs orption by 



the massive cool wind (e.g., Shore et al 



1996D from the 
ikolajewska 



RG component as suggested by Anupama & Mil 
(1999). 

Our ejected mass (1.9 x 10 _6 M Q ) is consistent with the 
observed ejected mass of (1.2 — 1.8)x10 _6 Mq estimated 
by Bohigas et al. (1989) except their distance of 1.6 — 2.0 
kpc to RS Oph. Using their equation (8) and V s = 200 
and V e — 2000 km s _1 , we obtain the cool wind mass of 
6 x 10~ 5 Af Q , indicating the wind mass loss rate from the 
RG of 3 x 10 _6 Mq yr -1 . It is reasonable that abo ut one 



twentieth of the cool wind accretes to the WD (e.g., Livic 



et al. 1986a) to power the TNR explosion. 



Hjellming et al. (1986) estimated the distance to RS 
Oph of = 1.6 kpc from H I absorption line mea- 
surements, using the hydrogen column density of Nn = 
(2.4 ± 0.6) x 10 21 cm~ 2 and the relation of N n /(T S d) = 
1.59 x 10 19 cm- 2 FT 1 kpc" 1 together with T s = 100 K. If a 
large part of this hydrogen column density stems not from 
the Galactic one but from the local one belonging to RS 
Oph, the distance is overestimated. For example, in the di- 
rection of the recurrent nova U Sco, Kahabka et al. (1999) 
reported a hydrogen column density of (3 — 4)xl0 21 cm~ 2 , 
which is much higher than the Galactic one of 1.4 x 10 21 
cm -2 . Here, U Sco is at least 2 kpc above the Galactic 
plane (b = 22°). If RS Oph is the same case as U Sco, the 
distance of 1.6 kpc is an upper limit and a much shorter 
distance is possible. 
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RS Oph 1985 Outburst (X=0.50, m =11.09) 




JD 2446000+ (day) 



Fig. 1. — Model light curves are plotted against time (JD 2,446,000+) together with the observational points of RS Oph outbursts. Filled 
circles indicate observational points with the previous outbursts being superposed (taken from Rosino 1987). Open circles correspond to the 
observational points of the 1985 outburst (taken from IAUC). Dashed lines denote the visual magnitude of the WD photosphere only, while 
solid lines represent the visual magnitude of the WD photosphere and the irradiated ACDK. The hydrogen content of the WD envelope is 
assumed to be X = 0.50. The apparent distance modulus is mo = 11.09. The mass of the WD is estimated to be 1.35 ± O.OIMq. 
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RS Oph 1985 Outburst 
1.35 M WD + 0.7 M RG m = 11.09 
T(RG)=3100 i=30° 7=0.4 X=0.50 

T(DK)=2000 a=0.01,0.10 0=0.01 v=9/8 7?=0.5 

WD + RG+AD with irradiation 
WD+RG with irradiation 

________ WD+RG with no irradiation 




wind stops 



end of steady H — burning 



1 00 



150 200 
JD 2446000+ (day) 
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Fig. 2. — Same as Fig. 1 but for the entire outburst phase. A dashcd-dotcd line is added and it denotes the visual magnitude of the WD 
photosphere and the irradiated RG companion. Thick and thin solid lines indicate the total visual magnitude of WD + RG + ACDK where 
the disk size is constant during the outburst, i.e., a = 0.01 and (3 = 0.01, and a = 0.01 and fi = 0.01, respectively. Two epochs are indicated 
by arrows with "wind stops" and "end of steady H-burning." The optically thick wind blows during the period from the very early phase of 
the outburst (HJD 2,446,092) to t ~ 67 days after maximum (HJD 2,446,159). The steady hydrogen shell-burning ends at t ~ 78 days after 
maximum (HJD 2,446,180). 
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RS Oph 1985 Outburst 

1.35 M WD + 0.7 M RG d=0.6 kpc (\=2.3) 
T(RG)=3100 i=30° 7=0.4 X=0. 35, 0.50, 0.70 
T(DK)=2000 a=0.01 0=0.01 u=9/8 7?=0.5 
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Fig. 3. — Same as Fig. 1 but for UV light curves. Three modeled light curves are calculated for X = 0.35, 0.50, and 0.70. Dashed lines 
denote total UV light from the WD photosphere, the ACDK, and the RG companion but dotted lines represent UV light only from the WD 
photosphere. Added large open circles are the observational points of ultraviolet (91lA — 3250A) observation by Snijdcr (1987). The distance 
to RS Oph is estimated to be 0.6 kpc if we assume an absorption of Ay = 2.3. These distance and absorption simultaneously satisfy both 
the V light curves (solid) and the UV light curves (dashed). 



